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ABSTRACT 

Issues regarding the use of optica! fiber thermometers to control heater settings in a microgravity 
vacuum furnace are addressed. It is desirable to use these probes in environments such as the International 
Space Station, because they can be operated without re-calibration for extended periods. However, the 
analysis presented in this paper shows that temperature readings obtained using optical fiber thermometers 
are corrupted due to emissions from the fiber when extended portions of the probe are exposed to elevated 
temperatures. 


INTRODUCTION 

The Quench Module Insert (QMI) is a platform for conducting solidification experiments in a 
microgravity environment. Current plans call for the installation of the QMI in the Microgravity Science 
Research Rank #1 on the International Space Station. The design of the QMI is similar to that of a 
Bridgman furnace and consists of a heater core, insulation jacket, instrumentation, coolant loop 
components, and a quench system. The heater core contains four heated zones and one water-cooled zone. 
This design produces the high thermal gradients required for directional solidification processing 
experiments. 

Precise thermal control is necessary to perform the candidate experiments for the QMI, and there is 
concern that thermocouples will drift due to exposure to cyclical thermal environments during the extended 
period of time it is planned to have the QMI in orbit. This concern led to the inclusion of optical fiber 
thermometers (OFT) in the preliminary design of the QMI. Preliminary designs call for the installation of 
Accufibers with sapphire blackbody sensors. Accufibers are a brand of OFT manufactured by LUXTRON 
Corporation. Product literature available from LUXTRON indicates that Accufibers have an accuracy of 
0.2% at 1000 C and a resolution of 0.01 C. In addition, Accufibers demonstrate excellent long-term 
stability and are immune to electromagnetic interference. In the proposed QMI design, the sapphire fiber is 
aligned with the axis of the QMI and an extended portion of the fiber is exposed to elevated temperatures. 
The analysis presented in this paper indicates that OFT readings will be corrupted under these conditions. 
Results obtained from a thermal model of the QMI are used to predict the temperature readings of an OFT 
probe. These predictions are consistent with readings obtained during the testing of a QMI prototype. 



OVERVIEW OF OPTICAL FIBER THERMOMETRY 


This section summarizes the operating principles of OFT (Dils 1983). An Accufiber sapphire 
blackbody sensor is illustrated in Fig. 1. The probe consists of a sapphire (AtOj) fiber whose sensing tip 
is given a metallic coating. The sensing end of the fiber is essentially an isothermal cavity, so the 
emission from this cavity will be approximately equal to the emission from a blackbody. The other end of 
the fiber is attached to the detection system. 



Figure 1 . Schematic diagram for an Accufiber sapphire blackbody sensor. The sensing tip (z=0) is coated 
with a thin metallic film to create a small isothermal cavity at a temperature of T 0 . The radiative flux 
emitted by the cavity, E^O), is approximately equal to the spectral emissive power of a blackbody, Eb X (' T 0 ). 


Modeling the sapphire fiber as a non-scattering medium, the radiative flux propagating along the 
fiber is governed by (Brewster 1992) 

^ = -K ai E, + K aX E bX (T(z)) (1) 

where E x is the spectral radiative flux (W/m 2 jim), is the spectral absorption coefficient (mm’ 1 ) and 
E bx (T(z)) is the spectral emissive power of a blackbody at a temperature of T(z). If the spectral absorption 
coefficient is independent of temperature, it is convenient to use the optical depth as the independent 
variable 




Eq. (1) then becomes 
^ + E x = E bl (T( tl )) 


( 2 ) 


( 3 ) 


Assuming that the sensing tip of the probe emits like a blackbody, the appropriate boundary 
condition for Eq. 3 is 

Ex(o) = E bX (r 0 ) (4) 

The solution to Eq. 3, subject to the boundary condition given by Eq. 4 is 


Ebx( T o) = E x (tx L )exp{t XL }- j E bX (T(t x ))exp{t x }dt x 


( 5 ) 


where t X L = K ax L. The spectral radiative flux measured by the Accufiber detection system, M x , is related to 
the spectral radiative flux at the end of the sapphire fiber. 

M x =E x (t XL )AXC (6) 


where AX is the width of the band pass filter used by the Accufiber probe and C is a correction factor to 
account for various losses in the detection system. Substitution of Eq. 6 into Eq. 5 gives 
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( 7 ) 


E b*(T 0 ) = ^^- J E b ,(T(tJ)exp{t x }dt x 

The integral on the right hand side of Eq. 7 represents the portion of the measured spectral 
radiative flux that is due to emission by the fiber. This integral can be neglected if sapphire is a poor 
emitter at the wavelengths of interest (K. x « 1) or if the fiber is at a low enough temperature to ensure that 
the emission by the fiber is negligible at the wavelengths of interest (T(%) « T 0 ). 

The spectral emissive power of a blackbody at a temperature of T is given by Planck’s equation. 



where the radiation constants are c, = 3.741 3x1 0 8 W|im 4 / m 2 and c 2 = I4388|imK (Brewster 1992). For 
the wavelengths used by the Accufiber probe (X { =0.80|um, X 2 = 0.95pm) and the temperature range of 

interest (400 - 1700 K), the exponential term is much larger than one. Therefore, Eq. 8 is well 
approximated using Wein’s limit. 



( 9 ) 


The temperature reading obtained by the Accufiber probe, T m , is calculated by neglecting the 
integral in Eq. 7. Substitution of Eq. 9 into Eq. 7 then gives 



X 5 M x q {xl 

c { AXC 


( 10 ) 


Using the measurements at two wavelengths and Eq. 10, the following ratio can be formed. 
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( 11 ) 


Assuming that AA,~AX 2 and that the loses in the Accufiber detection system do not depend on 
wavelength, T m is given by 



( 12 ) 


ERROR ANALYSIS 

Clearly, temperatures obtained using Eq. 12 will only be accurate when the approximations 
employed in the derivation are satisfied. In particular, emission by the fiber will change the value of the 
radiative flux measured by the Accufiber detection system and result in inaccurate measurements when large 
portions of the sapphire fiber are at elevated temperatures. Recall that the measured temperature is the 
temperature obtained when the integral in Eq. 7 is neglected. A relationship between the actual 
temperature, T 0 , and the measured temperature, T m , is obtained using Eq. 7 and Eq. 10. 
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( 13 ) 
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where 


f(*.T(t x )) = - J E bX (T(t x ))exp{t x ]dt 1 (14) 

c i o 

and T(tQ is the temperature profile along the sapphire fiber. 

In order to evaluate Eq. 13 and asses the errors due to emission by the fiber, it is necessary to 
know the temperature profile along the fiber and the absorption coefficient of AI2O3 at the wavelengths of 
interest. The following sections address these matters. 


TEMPERATURE PROFILE ALONG THE OPTICAL FIBER THERMOMETER 

The OFT probes are inserted in a boron nitride sleeve that surrounds the heated core of the QMI. 
The probes are coupled radiatively to the sleeve in which they are housed. There is also a conduction path 
along the fibers to a plate, which is used to support the probes. A SINDA model of the thermal 
environment of each fiber was developed to estimate the temperature profile along the fiber. A description 
of the geometry of the QMI and the SINDA thermal model is given in the Thermal Design Data Book 
(NASA-MSFC et al 1999) which was prepared for the Preliminary Design Review of the QMI. 

A prototype version of the QMI was tested at various heater settings. OFT readings were obtained 
for probes that are aligned with three of the QMI heating elements for the 6 cases listed in Table 1, The set 
points refer to the temperatures settings of each of the four heating elements, and these values are 10 to 40 C 
higher than the sleeve temperatures measured by the OFT probes. 


Table 1 . QMI Heater Set Points 


Case 

■E5E5EIIK3H 


HSSiSSmSH 


1 

600 

600 

600 

600 

2 

600 

600 

650 

650 

3 

900 

900 

900 

900 

4 

900 

900 

950 

950 

5 

1100 

1100 

1100 

1100 

6 

1200 

1150 

1100 

1100 


The SINDA model was run for each case listed in Table 1. Figures 2-4 show the estimated 
temperature profiles of the 3 probes. 



Figure 2. Estimated temperature profiles for Probe 1. 


TFAWS 99 


4 































1200 



Figure 3. Estimated temperature profiles for Probe 2. 



Figure 4. Estimated temperature profiles for Probe 3. 


OPTICAL PROPERTIES OF SAPPHIRE 

Brewster (1992) lists values for the real and imaginary parts of the refractive index a function cf 
wavelength. At 1 pirn, the imaginary part of the refractive index is k = 6 x 10' 8 . Assuming that the 
refractive index does not vary significantly with wavelength, the absorption coefficients at 0.80 and 0,95 
jim can be calculated. 


, = 47lk = 9.42x1 O' 4 mm' 1 
X, 

(15) 

_ 47tk ^ 94xio~’mm' 1 
112 X 7 

(16) 


These values are consistent with the data published by Gryvnak and Burch (1965) as shown in Fig. 3. 
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Figure 3. Spectral absorption coefficient of single crystal sapphire (AtOj) at elevated temperatures. From 
Gryvnak and Burch (1965). 


COMPARISON OF SENSING TIP TEMPERATURES WITH OFT READINGS 

Using the estimated temperature profiles shown in Figs. 2-4 and the spectral absorption 
coefficients given by Eq. 15 and Eq. 16, the integral in Eq. 14 was evaluated numerically. Estimates for 
T m were then calculated for each of the OFT temperature profiles. These values are compared with the 
measured OFT readings and the estimated temperature at the sensing tip of the OFT in Table 2. 


Table 2. Predicted and Measured OFT Readings 


Probe 

Case 

OFT Readings (C) 
Predicted Measured 

OFT Sensing Tip Temperature (C) 

1 

1 

598 

586 

594 ~ 1 

1 

2 

616 

591 

604 

l 

3 

900 

873 

890 

1 

4 

916 

879 

896 

1 

5 

1102 

1064 

1087 

I 

6 

1171 

1 132 

1167 

2 

1 

595 

612 

592 

2 

2 

621 

639 

612 

2 

3 

898 

910 

891 

2 

4 

920 

938 

903 

2 

5 

1100 

1 115 

1087 

2 

6 

1145 

1146 

I 142 

3 

1 

563 

576 

561 

! 3 

2 

618 

620 

616 

3 

3 

869 

857 

866 

3 

4 

924 

904 

920 | 

3 

5 

1075 

1052 

1071 

3 

6 

1085 

1059 

1080 


Comparison of the predicted OFT readings with the estimated sensing tip temperatures indicates 
that errors due to fiber emission increase as the length of fiber exposed to elevated temperatures increases. 
This results also show that the errors increase as the temperatures increase. 

The agreement between the measured and predicted OFT readings is somewhat imprecise. The 
measurements obtained from probe 1 are consistently lower than the predicted values. Measurements 
obtained using probe 2 are consistently higher than the predicted values. Since the temperatures in the 
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boron nitride sleeve vary considerably, these results indicate that the location of the probe’s sensing tips 
in the QMI prototype may differ from their location in the SINDA model. In addition, uncertainties 
regarding the thermal coupling between the OFT probes and their environment in the SINDA model make 
a precise interpretation of these results difficult. 


SUMMARY AND CONCLUSIONS 

The analysis presented in this paper demonstrates that an elevated temperature profile along the 
fiber will increase the temperature read by an OFT system due to radiative emission by the fiber. An 
expression that quantifies the errors due to fiber emission has been derived. The results presented in this 
paper show that the difference between the measured and tip temperatures decreases as the length of fiber 
exposed to elevated temperatures decreases. Also, the errors generally increase as temperatures increase. 

Predictions of the OFT readings based on a SINDA model of the probes thermal environment were 
compared with OFT readings obtained during the testing of a prototype of the QMI. The agreement 
between the predicted and measured values is not exact, but is consistent with uncertainties regarding the 
exact position of the sensing tips of the probes and the thermal coupling between the probes and their 
surroundings. Efforts to more accurately characterize the thermal environment of the OFT are currently 
being made. 
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